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ABSTRACT 



Geminga and B0656+14 are the closest pulsars with characteristic ages in the 
range of 100 kyr to 1 Myr. They both have spindown powers of the order 3xl0 34 
erg/s at present. The winds of these pulsars had most probably powered pulsar 
wind nebulae (PWNe) that broke up less than about 100 kyr after the birth of the 
pulsars. Assuming that leptonic particles accelerated by the pulsars were confined 
in the PWNe and were released into the interstellar medium (ISM) on breakup of 
the PWNe, we show that, depending on the pulsar parameters, both pulsars make 
a non-negligible contribution to the local cosmic ray (CR) positron spectrum, 
and they may be the main contributors above several GeV. The relatively small 
angular distance between Geminga and B0656+14 thus implies an anisotropy in 
the local CR positron flux at these energies. We calculate the contribution of 
these pulsars to the locally observed CR electron and positron spectra depending 
on the pulsar birth period and the magnitude of the local CR diffusion coefficient. 
We further give an estimate of the expected anisotropy in the local CR positron 
flux. Our calculations show that within the framework of our model, the local 
CR positron spectrum imposes constraints on pulsar parameters for Geminga 
and B0656+14, notably the pulsar period at birth, and also the local interstellar 
diffusion coefficient for CR leptons. 

Subject headings: pulsars: individual (Geminga, B0656+14), acceleration of par- 
ticles, cosmic rays, diffusion 



1. Introduction 



Geminga and B0656+14 are the close st pulsars with characteristic ages in the range of 
100 kyr to 1 Myr jManchester et alil2005T ). They both have spindown powers of the order 



3xl0 34 erg/s at present. The winds of these pulsars had most probably powered pulsar wind 
nebulae (PWNe) that broke up less than about 100 kyr after the birth of the pulsars. The 
reason for this statement is that we do not observe PWNe associated with pulsars older than 
100 kyr. 

Assuming that leptonic particles accelerated by the pulsars prior to breakup were con- 
fined in the PWNe and were released into the interstellar medium (ISM) on breakup of the 
PWNe (i.e. when the pressure of the ambient ISM exceeds the magnetic pressure B 2 / (8tt) 
in the PWN), we calculate the contribution of these particles to the locally observed cosmic 
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ray (CR) electron and positron spectra. We further calculate the expected anisotropy in the 
positron local interstellar spectrum (LIS) in the case of energy-dependent diffusion. 



Positrons from PWNe 



We discuss the acceleration o f particles by pulsars in the framework of the polar cap 
(PC) model (see e.g. the review of iBaring! 120041 ) . Given the relatively high surface magnetic 
fields of B s = 1.6el2G and B s = 4.7el2G for Geminga and B0656+14 respectively, a single 
primary electron released from the stellar surface will induce a cascade of electron-positron 
pairs in the magnetospheres of these pulsars. This amplification is modeled by introducing a 
multiplicity M' . The flux of primary electrons from the pulsar PC is given by the Goldreich- 
Julian current: 



I GJ w 2cp GJ A PC 



(1) 



with pgj the Goldreich- Julian charge density (jGoldreich &: Julian! Il969l ). and Ape ~ n Elp C 
the PC area. 

We assume that the electrons and positrons from the pulsar are re-accelerated at the 
pulsar shock, and model the particle spectrum by a p ower law with spectral index of 2, and 
with a maximum energy of ( Venter fc de Jagerll200( 



EL 



ecu 




(2) 



where k is the c ompression ratio at the sh ock, L S( j the spindown power, and a the magnetiza- 
tion parameter (IKennel fc Coronitilll984l ). This maximum energy stems from a condition on 
particle confinement: we require that the ratio between the Lamor radius and the radius r$ 
of the pulsar shock should be less than unity. We assume a maximum ratio e = r^/rs = 0.1 
and that for larger Lamor radii, the curvature of the shock results in particle losses. For 
the calculations presented in this paper, we assume k = 3. As observat ions suggest that a 
depe nds on the age of the pulsar (er = 0.003 is found for the C rab pulsar (jKennel fe Coroniti 
19841 ). which has an age of lkyr, but o = 0.1 for llkyr Vela (ISefako fc de Jagerl 120031 )). we 
thus assume 

/ t \ 3/2 

a(t) = .o^J , (3) 

with cr = 0.003. The fraction of the spindo wn power deposited in particles can be written 
in terms of the magnetization parameter a (IBuesching et al.l 120061 ): 

1 



?7part 



./ 



a 



(4) 
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Here, we int roduce a geometry factor / = VL acc / (47r), as a is supposed to be significantly less 
than unity (IKennel Coronitilll984l ). implying a large ^ part - Without /, from Eq. H]?7p ar t is 
of the order unity for a <C 1. This however is in contradiction with observations, indicating 
^ part < 1, as found e.g. for the Vela PWN JSefako fc de Jagerll2003h . We adopt / = 1/2. 



At any time the conditions 



M'L 



G.J 



Q'(E,t)dE 



Q'(E,t)EdE = r) p3It L sd (t) 



(5) 
(6) 



have to hold (jSefako fc de Jagerl 120031 ) . where Q'(E,t) = K'E~ 2 is the assumed particle 
spectrum at the pulsar wind shock. E m \ n is assumed to be similar to the inferred value of 
the Crab. A value of 100 MeV is adopted. 



For a non- decaying pulsar mag netic field (i.e. PP n 2 
L sd is given by iRees fc Gunnl (119741 ) 



PoPo 



n-2\ 



L sd {t) 
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sd,0 
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n + l 
' n—1 



the time-evolution of 



(7) 



with n = 3 representing a dipolar magnetic field, and r = Pq/ ((n—l)P ), P Q the birth period, 
and P the period's time-derivative at pulsar birth. (The subscript '0' denotes quantities at 
pulsar birth). These quantities are connected to the spindown luminosity at birth via 



J sd,0 



47T 2 /P n 



p3 



We assume that the particles are confined in the PWN for a time T, after which the PWN 
breaks up and releases them into the surrounding ISM. The time evolution of the particle 
spectrum Q(E,t) inside the PWN is described by 



d -^§^ - Q'(E,t) = A (B PWN (tfE 2 Q(E,t)) 



where we assumed a decaying magnetic field in the PWN 

1200 



B 



PWN 



(t) 



t/kjry 



(9) 



(10) 



This parametrization for B is justified as follows. After 10 to 20 kyr, the PWN field strength 
is already of the order of 5 /iG as observed by HESS from a number of PWN (see e.g. Ide Jager 
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20081 ) and after ~ 100 kyr we expect that the ISM pressure will randomize the PWN field 



structure, leading to relatively fast diffusive escape of charged particles from the nebula. 
Although the actual age for breakup is difficult to estimate, we assume to a first order a 
number of less than 100 kyr. GLAST obs ervations of a l imiting age for mature PWN may 
shed more light on this epoch of breakup (jde Jagerl 120081 ) . 

Eq. [9] can be solved using the Green's function formalism. The particle spectrum at 
time T is given by 

Q(E,T) = [ Q'(E ,t )E^E- 2 e(E - £ min ) 6 (£ max - E )dt , (11) 
Jo 



where O is the Heaviside step function and 

E 

E j t t0 B PW v(t>) 2 dt' + 1 



Our model predicts ~equal numbers of positrons and electrons to be accelerated by pulsars, 
thus Eq. [TT] also describes the source function for CR electrons. 



3. Propagation of CR positrons and local anisotropy 

The propagation of CR electrons and positrons in case of the diffusion coefficient k being 
spatially constant, is described by 

where N is the differential number density, S the source term, A the Lapace operator, and 
b the rate of energy losses. For a functional form of the diffusion coefficient k = koE a , with 
a = 3/5 and the energy losses b = boE 2 (i.e. synchrotron and inverse Compton losse s) one 
can find a Green's function solving Eq. (fT3"j) in the literature (IBerezinskii et al.lll990l ). It is 
given by 



G{f,fo,E,E ,t,to) = sft-to-^-iE^-E- 1 )^ 



exp 



(r-rp) 2 



6(ttA) l 



(14) 



with 



k (E a - 1 - E% ) , x 

A = 4-?A- '- . 15 

l-afco 1 ^ 
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Thus, the solution of Eq. ffl3|) is given by 



N 




G(r, r , E, E , t, t )Q(E )6 (t - U) 5 (r - fj) dE dt dr , 



(16) 



where Q(E ) is given by Eq. ffTTj) . fl and £, are the place and time of injection, repectively. 
We note, that the energy E of a particle at time t is linked to its energy E Q at injection by 



E = E/((t-t i )b E + l). 



(17) 



The anisotropy in the CR flux can be calculated in the context of diffusion as (IGinzburg fc Syrovatskii 
19641 ) 



5 



^max H~ -^min 



3£;|ViV| 



18) 



where VA" denotes the gradient of AT. The expected anisotropy in the positron LIS was 
calculated assuming the contribution of a nearby source, as given by Eq. ([TBI , on top of an 
i sotrop ic background. For the background we assumed a power law fit given by lBarwick et al. 



(119981 ) . The calculated anisotropies are given in the right panels of Figs. [T] and [2] (thin lines). 



To get an estimate of the maximum expected anisotropy, we also calculate the anisotropy 
assuming that the whole CR positron flux originates from a po int source (th i ck lin es in the 
right panels of Figs. [1] and [2]). For energy- independent diffusion, iMao fc Shenl (119721 ) derived 
the simple relation 



5= *A 
2c U 



Allowing for energy- dependent diffusion, we get, inserting Eq. ( Tl6|) into Eq. ( fl8l) . 



5 = —r t b (a - 1) E a (E^ 1 - E^ 1 ) 1 
— c 



(19) 



(20) 



As in the case of Eq. (Tl9|) . Eq. (TSOj) also does not depend on the magnitude of the diffusion 
coefficient, only on the distance to the source, and via Eq. (fTTj) . on the time since the 
injection of the particles into the ISM. In the limit of E — ► 0, Eq. ( f20l) reduces to Eq. ( fl9l) . 



We calculated t he contribution from Geminga an d B0656+14 to the pos itron LIS for 
distances of 157 pc (jCaraveo et al.l Il996l ) and 290 pc ([Manchester et al.l 120051 ) respectively 
(assuming T = 20kyr and T = 60kyr), in addition to the expected anisotropies in the 
positron LIS. The results for birth periods of 40 ms and 60 ms are plotted in t he left panels 



of Fig s. [T]and [21 where we compare our calculations with the measurements from lBoezio et al. 
(hoOOf ) and lDuVernois et al.l (bpOlh . 
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Fig. 1. — Left panel: Contribution of Geminga to the positron LIS for fc = 0.1 kpc 2 Myr -1 
and P = 40ms, T = 20kyr (long dashed line), P — 40ms, T = 60kyr (dot-dashed 
line), and Pq = 60ms, T = 20kyr (dashed) on top of an isotropic background (solid 
line). The thin lines mark the combined spectra (pulsar contribution plus background), 
whereas the t hick l ines give the cont ribution of the pu l sar al one. Also shown are data from 
Boezio et al.l (120001 ) (diamonds) and iDuVernois et al.l (120011 ) (triangles). Right panel: the 
expected local anisotropy in case only Geminga contributes to the LIS (thick lines), and in 
case Ge minga contributes on top of an isotropic background positron flux (thin lines) as 
given by iBarwick et al.l (119981 ) (solid line in left panel). The line styles correspond to the 
cases as given for the left panel. The thick dashed and long dashed lines coincide. 



Conclusions 



We have shown that one can expect a non-negligible CR positron component in the 
LIS from nea r by pu lsars that may become dominant above several GeV, in agreement with 
Atoyan et al.l (119951 ) who showed that the high-energy positron LIS may be explained by a 
young, nearby source. In the context of our model, we are able to constrain the permissible 
pulsar birth period P , depending on the magnitude of the interstellar diffusion coefficient. 
For the two nearest pulsars with characteristic ages in the range le5yr to le6yr, Geminga 
and B0656+14, we show that in particular for B0656+14 one can expect, in the absence of 
a background flux, an anisotropy in the positron LIS of up to almost 3%, significantly larger 
than the expected value of ~0.25% for Geminga. As shown in Figs. [I] and [2], the observed 
anisotropy also gives an estimate of the contribution of the pulsar to the positron LIS. On 
the other hand, a measured anisotropy larger than the ~3% we obtained for B0656+14 
would be, from Eq. (1T91 and Eq. (1211 . an indication of the existence of an even younger, 
nearby source, e.g. a longer lifetime T of the PWN. We also note that the predicted flux 
from PSR B0656+14 appears to overpredict the observed flux above 10 GeV. This implies 
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Fig. 2.— Same as Fig. Q] but for B0656+14. 



more severe constraints on the pulsar output, whereas Geminga's parameters are not that 
severely constrained by CR positron observations. 

We remark that Galactic CR, including electrons and positrons, are subjected to solar 
modulation at energies below «10 GeV. The encounter of these particles with the solar wind 
and imbedded magnetic field causes a heliospheric anisotropy that is primarily determined 
by the combined modulation effects of convection, diffusion, and drifts - all solar cycle 
dependent. Drifts will cause this anisotropy to have a 22-year cycle. CR electrons and 
positrons at 1 GeV at Earth may therefore exhibit a heliospher ic anisotropy of up to a few 
percent, assuming that they enter the heliosphere isotropically (IPotgieter fc Langnerl 12004 ). 
It will be an interesting exercise to determine how this anisotropy will change if the LIS is 
anisotropic. However, the anisotropy that we pre dict here is the lar gest above 10 GeV, an 
energy range at which only the PAMELA mission (IBoezio et al.ll2004j ) may be able to gather 
sufficient statistics to find an anisotropy of the predicted magnitude. 
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